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Abstract 

It is widely believed that dark matter exists within galaxies and clusters of 
galaxies. Under the assumption that this dark matter is composed of the lightest, 
stable supersymmetric particle, assumed to be the neutralino, the feasibility of its 
indirect detection via observations of a diffuse gamma-ray signal due to neutralino 
annihilation within M31 is examined. 
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1 Introduction 



The existence of cosmic dark matter (DM) is required by a multitude of observations 
such as excessive peculiar velocities of galaxies within clusters of galaxies or gravitational 
arcs. Furthermore, both Big Bang nucleosynthesis (which predicts a baryonic relic density 
Qb <S Cltat — 1) an d plausible scenario for large-scale structure formation, strongly suggest 
a substantial non-baryonic/cold DM component in the Universe. It so happens that 
super symmetric extensions of the standard model of particle physics provide a natural 
candidate for such a DM in the form of a stable uncharged Majorana fermion (Neutralino). 
Hereafter, we briefly report on the potential of the high-energy 7-rays (E 1 > 50 GeV) 
ground based detector CELESTE (de Naurois et ai. '02), to detect neutralinos indirectly 
through their annihilation in the halo of M31. (For the detailed study see Falvard et ai. 
'02) 



2 Neutralino halo around M31 — Neutralino anni- 
hilation 

The late-type Sb spiral galaxy M31 lying at a distance of 700 kpc has a visible part 
consisting mostly of a bulge and a disk. We have reconsidered the two mass components 
fit to the rotation curve of M31, performed by Braun ('91). Taking the following mass-to- 
light ratios T bu ige = 6.5 ± 0.4 T B and T disk = 6.4 ± 0.4 T Bi0 (where T B is the mass- 
to-light ratio for the Sun) Braun concluded that no dark halo is necessary to account for 
the velocity field. However, this conclusion relies on a large value of Tdisk which disagrees 
with estimates based on the blue color of the disk and on synthetic spectra of young stellar 
populations which it contains (Guiderdoni '87). The mass-to-light ratio Tdisk of a purely 
stellar component should actually not exceed ~ 3.8 Tb i0 . In addition, a disk as massive 
as that proposed by Braun should generally be unstable. Therefore, we have assumed 
the presence of an additional mass component in terms of a spherical halo whose mass 
density profile is generically given by 

r V f Tq + a c 



where 7, a and e define the various profiles. While a fit to the observations constrains 
weakly the ratio T d isk/T bulge, it leads to more stringent constraints on the structure of the 
neutralino halo, actually favouring a NFW profile 7=l,a = l,e = 2 (Navarro, Frenk, & 
White '96). Typical results are illustrated in Fig. 1. 

If neutralinos are indeed a substantial component of the DM around M31, one expects 
them to annihilate into standard model particles, eventually producing energetic photons. 
The corresponding photon flux at Earth - per unit of time and surface - may be expressed 
as 

I 1= ^h- = / / JdidO = i-<^E (2) 

dt dS 47r m 2 x Jfov J ios x 47r m? x 
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where m x is the neutralino mass and (av) iV 7 denotes the thermally averaged annihilation 
rate yielding ]V 7 gamma-rays in the final state. Eq. @ encapsulates all the particle physics 
features in (av) iV 7 /m x , while the astrophysical modelling is contained in the line of sight 
integral S. In Table 1 of Fig. 3 we illustrate some typical fluxes taking a nominal neutralino 
mass, m x = 500 GeV, and annihilation cross section, (av) iV 7 = 10~ 25 cm 3 s' 1 . 

3 Super symmetric model predictions and resulting 
fluxes 

In the present section we consider more specific particle physics model predictions of the 
7-ray fluxes. We will focus mainly on the minimal supergravity scenario (mSUGRA), 
(Barbieri et al. '82, Chamseddine et al. '82, Hall et al. '83). Making the usual sim- 
plifying assumption of common universal values at some grand unified theory (GUT) 
energy scale M GUT ~ 2 x 10 16 GeV, i.e. m sca i ars (M GUT ) = m , M gauginos (M GUT ) = 
mi/2, Atriiinear(MGUT) = A , and requiring various physical consistencies (e.g. electroweak 
symmetry breaking at low scale, electrically neutral and stable DM particles, present ex- 
perimental limits from particle colliders, etc..) one can make specific predictions for the 
7 flux in terms of these model parameters. Assuming the neutralinos should account for a 
large fraction of the DM at cosmological scales as well, one also requires the corresponding 
relic density to be in the right hallmark, 0.025 < Q x h 2 < 0.3. In Fig. 2 we present a 
scan over the basic parameters of mSUGRA for the correlation between the integrated 
7-flux and the neutralino mass and relic density. All results have been obtained using 
an interfaced version of the two public codes DarkSUSY and SUSPECT. Convoluting the 
annihilation spectra with the CELESTE acceptance, we show in Fig. 3 the expected rate of 
7/min from a maximal M31 smooth halo (see Fig.l). This puts the signal from such halos 
clearly beyond a reasonable CELESTE reach which is roughly one tenth of the CRAB 
nebula (4.9 7/min). Nonetheless, a significant departure from a smooth halo is not unre- 
alistic due to possible dumpiness and/or black hole accretion effects. This can lead to a 
global enhancement factor which may reach up to two orders of magnitude, as illustrated 
in Table 2., making the survey of M31 with CELESTE worth being undertaken. 

4 Conclusion 

We conclude that under favourable conditions such as rapid accretion of the neutralinos on 
the central black hole in M31 and/or excessive halo dumpiness, a neutralino annihilation 
7-ray signal may be seen by the ongoing observations of M31 with CELESTE. 
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Figure 1: A 7 = 1 neutralino halo is added to the bulge and to the disk of M31. Left: 
an intermediate case with T bu i ge = 4.2 Tb,q and T disk = 4.2 T B . Right: a maximal 
halo with T bulge = 3.5 T B , T disk = 2.5 T B . The global solid rotation curve is in good 
agreement with the data of Braun. 
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Figure 2: Left: the integrated 7 flux from M31 as a function of m x for i? 7 > 30 GeV. 
Each point corresponds to a model in our "wild scan". Three different ranges of gaugino 
fraction are considered. Right: the integrated 7 flux from M31 as a function of Q x h 2 . A 
renormalization procedure for Q x h 2 < 0.1 has been applied on flux values (see Falvard A. 
et al. '02). 



T bulge T** E„(3.E 


kpc) E I ,(28kp 


i / (3.S kpc) 


J, (28 kpc) 


6.5 6.4 











4.2 4.2 1 


1.2 


3.2 X 10-" 


3.8 X 10-" 


3.5 2.5 3 


3.7 


10.2 X 10-" 


11.8 X 10-" 


Table 1: 


Chimp distribution 


Dynamics 


SBH environment E S?" 


compact 


no evaporation 


SIS spike 


165 100 


extended 


no evaporation 


SIS spike 


140 90 


compact 


no evaporation 


p oc r~ 


80 20 


extended 


no evaporation 


p oc 


55 6 


compact 


evaporation 


p cx r -1 * 8 


25 6 


extended 


evaporation 


/i k r- 1 - 8 


40 5 




100 200 300 400 500 600 700 
Neutralino Mass (GeV) 



Figure 3: Table 1: Three different models for M31 are featured. The first corresponds to 
the case where no halo is needed. S 19 (i?) is in units of 10 19 GeV 2 cm~ 5 and the flux I 7 (R) 
is for a circular region encompassing the inner 3.5 kpc (corresponding to the CELESTE 
f.o.v. of 10 mrad) and 28 kpc, respectively. Table 2: Impact of astrophysical parameters 
such as dumpiness of the M31 halo and supermassive black hole (SBH) in its centre, 
on flux predictions, smooth and clump contributions added. (S 10 corresponds to the 
CELESTE f.o.v.). Right: Number of 7/mm as expected with CELESTE from M31. 
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